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We expose the relation between the properties of the three-body continuum states and their 
two-body subsystems. These properties refer to their bound and virtual states and resonances, 
all defined as poles of the 5'-matrix. For one infinitely heavy core and two non-interacting light 
particles, the complex energies of the three-body poles are the sum of the two two-body complex 
pole-energies. These generic relations are modified by center-of-mass effects which alone can produce 
a Borromean system. We show how the three-body states evolve in ®He, ®Li, and ®Be when the 
nucleon-nucleon interaction is continuously switched on. The schematic model is able to reproduce 
the main properties in their spectra. Realistic calculations for these nuclei are shown in detail for 
comparison. The implications of a core with non-zero spin are investigated and illustrated for ^’^Ne 
(^®0-|-p-|-p). Dimensionless units allow predictions for systems of different scales. 

PACS numbers: 21.45.-l-v, 31.15.Ja, 25.70.Ef ll.80.Jy 


I. INTRODUCTION 

The three-body problem has a long history from 
macroscopic celestial classical mechanics, e.g. sun-earth- 
moon Q to microscopic quantum mechanics, e.g. the 
helium atom 0, three nucleons Q or three quarks 0. 
The modern treatment was boosted by the formulation 
of the Faddeev equations originally aimed at scat¬ 
tering problems, but also successfully applied for bound 
states |a,l3- The more recent interest in bound state halo 
structures and Borromean systems are by now fairly well 
understood in terms of the basic two-body interactions 
13. The success is at least indisputable for three-body 
systems with only one or a few bound states. 

On the other hand, the corresponding three-body prop¬ 
erties for positive energies (energies above breakup) are 
much less established although studied thoroughly for 
both short and long range interactions i i. Con¬ 
tributions from both short and long-range interactions 
make computations numerically difficult. The three- 
bod y C oulomb problem itself is still considered unsolved 
[ifl fill fl3 and three-body resonances for stroi^y in¬ 
teracting particles are still debated Q, 0, 0, . 

This is unfortunate, since the continuum structure often 
form the basis in descriptions of the dynamic behavior of 
a given system. 

Crucial properties of the continuum are revealed by 
information about the resonances and virtual states de¬ 
fined as poles of the ^-matrix. Together with the discrete 
set of bound negative-energy states we also have a dis¬ 
crete set of unbound complex-energy states with their 
corresponding wave functions. For completeness also the 
continuous non-singular background states are needed, 
but the singularities of the scattering matrix are very of¬ 
ten decisive. Important examples using different methods 
within few-body physics are astrophysical reaction rates 
adiabatic reaction processes arising at low energies 


or for large impact parameters m, three-body decays 
[^ . three-body resonances for Faddeev operators 
for nucleons [ 23 |, for electrons and positrons and 
four-body nuclear continuum states |24l| . This fist could 
be extended. 

Continuum structure is in general more difficult than 
bound state problems although various methods have 
been designed to overcome the technical problems at least 
for resonances, see for example n m m m m 123 
An understanding of the generic origin of ^'-matrix poles 
would be tremendously helpful especially if recognizable 
traces of a well-structured origin are left in the realistic 
spectra. No doubt this would allow easier interpretation 
of complicated numerical results, allowed design of bet¬ 
ter methods, and indicate which effects to look for in 
different contexts. 

The purpose of this work is to investigate how the 
two-body interactions determine the three-body contin¬ 
uum structure. We focus on three-body continuum states 
where none of the two-body subsystems is bound (Bor¬ 
romean systems). When two-body bound states are pos¬ 
sible, different three-body structures can appear, as for 
instance, unbound three-body states with negative en¬ 
ergy. However, these systems can very often be under¬ 
stood as two-body systems made by a two-body bound 
state and the remaining third particle in the continuum. 
With the help of the hyperspheric adiabatic expansion 
method we determine which three-body bound or vir¬ 
tual states and resonances result from given sets of cor¬ 
responding two-body properties. Taking a simple sys¬ 
tem as starting point, it is then possible to observe how 
the three-body states evolve in the complex energy plane 
when more and more realistic features are incorporated 
into the calculation. 

In section II we give details about the complex 
scaled adiabatic expansion methods. In section III the 
schematic system of an infinitely heavy core and two mu- 
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tually non-interacting light particles is described analyt¬ 
ically, while in section IV the general properties of the 
system are described after numerical studies of specific 
systems. In section V we relate the spectra from the 
schematic model with those obtained in realistic calcula¬ 
tions for systems with zero (®He, ®Li, ®Be) and non-zero 
core-spin (^^Ne). We finish in section VI with some qual¬ 
itative estimates for other systems and in section VII we 
give a short summary and the conclusions. 


II. THE COMPLEX SCALED HYPERSPHERIC 
ADIABATIC EXPANSION METHOD 


To describe a three-body system we use the standard 
coordinates: 
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where mi, m^, and m^ are the masses of the three par¬ 
ticles and m is an arbitrary normalization mass. From 
the Jacobi coordinates we define the hyperspheric coordi¬ 
nates {p, fij}, where {fl*} = {ai, Oa,,, p=\/xl + y'f, 
ai=&vcta,n{xi/Ui)^ and and give the directions of 
Xi and respectively. 

The three-body wave function is written as a sum of 
three components tl;^'‘\xi, y^), each corresponding to one 
of the three possible sets of Jacobi coordinates 0. These 
three components satisfy the three Faddeev equations 

(r-F;)V-«(a;„y,) + V,fe(crO (2) 

X (^i}^'^\x,,y^)+i}^^\x^,yj) + 'ij;^'^\xk,yk)) =0, 

where T is the kinetic energy operator, Vjk(xi) is the two- 
body interaction between particles j and k, and E is the 
three-body energy. Here {i,j,k) is a cyclic permutation 
of (1,2,3). 

We employ the hyperspheric coordinates to expand 
each component {xi,yj) in terms of a complete set 
of angular functions 

^ ^(3) 

^ n 


where the additional index no labels different solutions 
we later on want to consider. Usually the expansion 
converges rather fast, and only a few terms (typically no 
more than three) are needed. 

By rewriting Eq. in terms of the hyperspheric coor¬ 
dinates, and inserting the expansions in Eq. © we sepa¬ 


rate the Faddeev equations into angular and radial parts: 



n' ^ ^ ' 


where is a hyperangular operator which together with 
expressions for the functions (p) and Qnn' (p) can be 

found in 0. The complete set of angular functions 
in the expansion © is precisely the eigenfunctions of the 
angular part of the Faddeev equations. The index n la¬ 
bels the the corresponding eigenvalue A„, which enters in 
the coupled set of radial equations © as effective poten¬ 
tials. The index uq is related to the boundary condition 
for continuum wave functions, which for short-range in¬ 
teractions and no two-body bound states has the asymp¬ 
totics given by 0 : 


where no then labels the incoming channel, the S'-matrix 
is a factor depending on the complex momentum k 
related to the complex energy E by K=-^2mE/h?, and 

and are Hankel functions of first and second 
kind. The hypermomentum K is given by the asymp¬ 
totic behaviour of the angular eigenvalue A„(p) that ap¬ 
proaches K{K + 4) 0 . 

For three-body bound states must fall off expo¬ 
nentially at large distances. Then only one Hankel func¬ 
tion is present and the S'-matrix S"“ has a pole for the 
imaginary momentum k {k = f|K|) and the energy E 
is negative. Poles of the S-matrix in the lower half of 
the complex momentum plane correspond to three-body 
resonances, and their asymptotic radial wave function is 
given only by the in Eq.®. Asymptotically 

the Hankel function FJA+2 vanishes exponentially like 

g-Elp^ while H^x +2 grows like This means that the 
radial coefficients of the continuum wave functions are 
dominated by the PEW't. Therefore, it is very dif¬ 

ficult to distinguish the resonance wave function, where 
only the form is present in the asymptotics, from 

an ordinary continuum wave function. 

This problem is solved by applying the complex scaling 
method |2^ , where the Jacobi coordinates 

Xi and yi are rotated into the complex plane by an arbi¬ 
trary angle 6 {xi yi j/iC*®). This means that 

only the hyperradius p is rotated {p pe*®), while the 
five hyperangles {12^} remain unchanged. After this ro¬ 
tation the radial wave functions of the resonances behave 
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asymptotically like 


_j g-|K|psin {9-9ft) ^i(\K\pcos {9-9ft.)-KTT/2+3Tr/4:) 

where On is the argument of the complex momentum k 
[k = |K|e“*®®). From o we observe that when 9 > Or 
the radial wave function falls off exponentially, exactly as 
a bound state. Continuum wave functions are dominated 
in the asymptotics by the term, that 

diverges exponentially when 9 > Or. Thus, after com¬ 
plex scaling, resonances can be easily distinguished from 
ordinary continuum states, and furthermore, the numer¬ 
ical technique used to compute bound states can be used 
for resonances. In particular, after solving the complex 
scaled equations and 0 with the boundary condition 
o three-body resonances and bound states are simulta¬ 
neously obtained. Solving Eq. with a box boundary 
condition {f{pmax)=^, Pmax being a large value of the hy¬ 
perradius) discretizes the continuum spectrum, and the 
continuum states are rotated by an angle 9 in the mo¬ 
mentum plane and 29 in the energy plane . The box 
boundary condition is often enough to obtain accurate 
bound state and resonance solutions. 


III. TWO INDEPENDENT SUBSYSTEMS 

We analyze in this section the schematic model of an 
infinitely heavy core and two light particles. These two 
mutually non-interacting light particles interact with the 
core via short-range interactions. We shall first assume 
that these three particles have zero spin and then later 
generalize to non-zero particle spins. 


A. Particles without spin 


From the Jacobi coordinates o we can construct the 
corresponding conjugated momenta: 
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In precise analogy to the hyperradial coordinates, we 
now introduce the hypermomentum coordinates {a, 

defined as {IIkJ = {a„,, k = ^k^. + k ^., 

= arctan(fca:^/fcyj. The three-body momentum k 
defined earlier is independent of the Jacobi coordinate 
system labeled i. The directions of fej,. and ky. are de¬ 
scribed by flk^. and , respectively. 

When one of the three spinless particles has infinite 
mass (let us take mfc=oo), Eqs.flJ and © simplify, in 
the three-body center-of-mass system where rk=0 and 
Pi+Pj+Pk=0, to the form in Fig^ Therefore, in two of 



FIG. 1: Sketch of the schematic three-body system. Particles 
i and j do not interact with each other. Particle k is infinitely 
heavy. One set of Jacobi coordinates and their corresponding 
conjugate momenta are given in the lower part. 


the three possible sets of Jacobi coordinates, the coordi¬ 
nates are simply proportional to the distances between 
the infinitely heavy core and each of the two light par¬ 
ticles. If the two-body interactions depend only on the 
distance between particles, the three-body Hamiltonian 
H can be written as H = Hki + Hkj (see FigCJ where 
Hki and Hkj are the two-body Hamiltonians describing 
the corresponding two-body systems. From this separa¬ 
bility, the three-body bound state energies are obviously 
given by the sum of the corresponding two-body bound 
state energies. 

For three-body resonances and virtual states the same 
happens. Three-body resonances show up at energies 
equal to the sum of the two-body resonance energies 
of the ki and kj subsystems. From the separability of 
the Hamiltonian this result may be accepted as obvious. 
However, the only trivial deduction using the separabil¬ 
ity is that the three-body energy is given by the sum 
of the two two-body energies. This does not necessar¬ 
ily imply that a three-body energy equal to the sum of 
two two-body resonance energies must correspond to a 
three-body resonance. Bound states have discrete en¬ 
ergies, and the split of a given three-body bound state 
energy into two two-body energies can be made only in 
specific ways, since the two-body energies can have only 
the specific values corresponding to the discrete two-body 
bound state spectrum. However, resonances are contin¬ 
uum states, and then a given three-body energy can be 
obtained by infinitely many pairs of two-body energies. 
Therefore, if a given three-body energy matches with the 
sum of two two-body resonance energies, there are also 
infinitely many other pairs of continuum (non-resonant) 
two-body energies whose sum gives the same three-body 
energy. Thus, it is not a trivial conclusion that only 
the particular three-body state of matching energy cor- 
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responds to a three-body resonance defined as a pole in 
the three-body S'-matrix. 

Let us then investigate this more closely. From the 
separability of the Hamiltonian it follows that the three- 
body wave function is the product of the wave functions 
describing the ki and kj two-body subsystems. For un¬ 
bound two-body systems the two-body wave function can 
be expanded in partial waves, and the three-body wave 
function takes then the form: 

4'(+)(a;,y, fey) = (9) 


the spherical harmonics and the Jacobi polynomials 0. 
Due to the delta functions, the three-body hyperradial 
functions in Eq. m are block-wise diagonal and each 
block matrix can, for given {lx,iy,L), be labeled by the 
indexes KK'. Note also that the wave function does not 
depend on the total orbital angular momentum L. All 
L-values allowed by coupling are degenerate. 

In the limit of zero interactions between the core and 
the light particles the two-body radial wave functions 
ui{x, k)/xk become the spherical Bessel function ji{xk), 
and Eq. 03 becomes: 



J Kl^lyL 


(kp) = 


'5LL'5KK'^ 


K 


Jk+ 2 {kp) 

[kpY ’ 

(13) 


where Jk +2 is a Bessel function of the first kind. Insert¬ 
ing (ini) into 03 we then obtain the three-body wave 
function 


where {x, y) refer to the Jacobi system in Fig^and the 
momenta of the corresponding subsystems are and 
fey. The spherical harmonics and the two-body radial 
functions are denoted Yimi and ui. 

Following M , also the continuum wave function for a 
three-body system of spinless particles can be written as: 


4'(+)(a;,y,fe,j,fey) = 

/ . ^ KLM 
KLMl^ly 


JKl^lyL 

K'L'M'l'i' 


( 10 ) 

{^p)Yk'um' (^) 


where {x,y), (fe^,, ky), and could correspond to any 
of the three possible sets of Jacobi systems. We omitted 

the index i. The hyperradial functions, \ are 

solutions to Eq. m where we now specified the indices 
n and Uq and explicitly included the factor k in the ar¬ 
gument. The hyperspheric harmonics can for 

instance be found in 0. Using the orthogonality of the 
hyperspheric harmonics, we immediately have 


AK't'jyL') 

JKixiyL 


(kp) = 


(11) 
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and by inserting Eq. into m we get 


d'(+^(a;,y,fe,r,fcy) = 
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and the plane waves are recovered from the hyperradial 
partial wave expansion. 

The large-distance asymptotic behaviour of the three- 
body radial wave functions, Eq. 03, can be obtained 
from the asymptotics of the two-body wave functions, 
that for the case in which none of the two-body core¬ 
particle subsystems is bound has the form: 


^ ^ (hf\xk) + Si{k)h[^\xk)^ , (15) 

(i) 

where fe) ^ are the corresponding two-body Hankel func¬ 
tions of first and second kind, and Si{k) is the two-body 
S'-matrix. Inserting Eq. 03 into 03 and using the ex¬ 
pressions given in the appendix, we obtain 


AK't'j'yL') 2 ^ 

fKtxtyL i'^P) — ^exi'x:^iyt'y^LL' I NNj^, 


Jo 

X f da{s\naY'‘^'^{cosaY'>^'^pl^Y^^\cos2a) 
Jo 


Ui^{x,kx) ue^{y,ky) 
xkx yky 


( 12 ) 


where is the normalization constant appearing 

when expressing the hyperspheric harmonics in terms of 


where we omitted the unimportant block indexes 
{ix,^y,L). Then Ea. (nT)ll is identical to Eq. after divi¬ 
sion by the phase factor 1 in Eq. ©. The three-body 
S-matrix takes the form: 


(^) 


2JY^‘^ da^<^>Yn'''\o^i^)six:{kx)siy{ky) 
da^^Yn'''\o:^)isixikx) + Siy{ky)) 


(17) 


where 


= (sinQ;„)^^"+2(cosa«,)2^’'+2 


(18) 


xPn ^ "'(cos2a^)P^, "'(cos2aK), 
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is a smooth function of given in terms of the Jacobi 
functions ^ ^ with K = 2n + and K' = 

2ti' P ix P ^y- 

The integrands in Eq. m contain the two-body S- 
matrices depending on corresponding momenta. The in¬ 
tegration variable maintains the total energy of the two 
subsystems while varying the ratio of their momenta. As¬ 
sume that one pole of one of the subsystems is a first 
order pole which does not coincide with any pole of the 
other subsystem, then the integrations in Eq. 03 across 
the pole give a smooth function without trace of the pole. 
However, if two poles from different subsystems exist si¬ 
multaneously the product of two-body S'-matrices in the 
numerator of Eq. 03 results in a pole term of second 
order, which after integration across the pole, leaves a 
first order pole in the three-body 5'-matrix. The energy 
of this three-body pole corresponds then precisely to the 
sum of the energies of the two two-body poles. 

If one pole is of higher order, and not coinciding with 
a pole from the other subsystem, it would survive the 
integration in precisely in the same way in numerator 
and denominator. Thus no trace is left in the three-body 
5'-matrix. A three-body pole only arises when two two- 
body poles belonging to different subsystems exist simul¬ 
taneously. Combining any order of coinciding poles from 
the two subsystems is then seen to produce a three-body 
pole of an order equal to the smallest order of the two- 
body poles involved. This proves that the three-body 
poles appear if and only if the energy is equal to the sum 
of energies of two two-body poles from different subsys¬ 
tems. The smallest order of the two coinciding poles 
reappears in the three-body pole. 

In case of having a two-body bound state, the corre¬ 
sponding two-body asymptotics is not given by Eq. m 
but by a falling off exponential e where k is deter¬ 
mined by the two-body binding energy. The three-body 
properties are then directly dictated by the two-body S- 
matrix describing the remaining core-particle subsystem. 
If this two-body S'-matrix has a pole corresponding to 
a two-body bound state, the three-body system is ob¬ 
viously bound with a binding energy equal to the sum 
of the two two-body binding energies. If the two-body 
S'-matrix has a pole corresponding to a two-body reso¬ 
nance (or virtual state), the three-body system has then 
a resonance (or virtual state) that actually corresponds 
to a resonance (or virtual state) of the two-body system 
made by the bound two-body state and the remaining 
third particle. 

At this point it is also important to mention that we are 
identifying poles of the S-matrix in the fourth quadrant 
of the energy plane and resonances. However, strictly 
speaking, for a pole to be considered a resonance it is 
required that its width is clearly smaller than the energy 
separation between resonances. In this way resonances 
are the poles of the S'-matrix close enough to the real 
energy axis. It can certainly happen that when two com¬ 
plex energies corresponding to two singularities in the 
S-matrix are added, the final energy can be close to the 


negative imaginary axis, or even in the third quadrant 
of the energy plane. These singularities should not be 
considered as resonances in the sense given above. 

The formalism leading to the expression for the three- 
body S-matrix is applicable for a very different system, 
i.e. three particles of finite mass with only one non¬ 
vanishing two-body interaction. If we choose the Ja¬ 
cobi coordinate system where x is related to the non-zero 
two-body interaction the wave function is again given by 
Eq.@ with a Bessel function instead of ut^/{yky). All 
derivations remain unchanged and we arrive at the S- 
matrix in Eq. 113 . The difference is that now only the 
two-body S-matrix related to the a;-degree of freedom has 
poles. Therefore the three-body system has no S-matrix 
poles. 

This system effectively also consists of two independent 
two-body systems, but now reached in a completely dif¬ 
ferent limit with one non-zero interaction and arbitrary 
masses. The result then illustrates other aspects of the 
continuum properties for three-body systems. In both 
cases the conclusion is that one subsystem alone cannot 
produce three-body resonances or virtual states. At least 
two different subsystems must collaborate and simulta¬ 
neously contribute with coinciding poles. It then seems 
inevitable that, if all the three two-body subsystems only 
have poles, then fully realistic three-body systems must 
also have only poles. 


B. Particles with spin 

The results in the previous subsection remain valid for 
particles of non-zero spin if the two-body interactions are 
spin-independent, because then the spin part of the wave 
function is decoupled from the coordinate part. In fact, 
even in the case of identical particles the spin part can 
always be used to establish the proper (anti)symmetry of 
the wave function without changing the coordinate part. 
The only exception is two identical bosons with zero spin 
where the spin part of the wave function always is sym¬ 
metric under exchange of the two bosons. Therefore odd 
values of the relative orbital angular momentum between 
the bosons are strictly forbidden. This fact leads to im¬ 
portant differences compared to non-identical particles 

When the two-body interactions are spin-dependent, 
the separability of the three-body Hamiltonian for the 
system shown in Figni can disappear. We denote the 
spins of particles i, j, and k by Sj, s^-, and s^, respec¬ 
tively. Let us for simplicity assume that the particle-core 
interactions of the system in FiglD contain a central plus 
a spin-spin term, i.e., Vki = ■ Sfe, where 

ji = f-ki + Si and tki is the relative orbital angular mo¬ 
mentum between the core and particle i (and similarly for 
the Vkj interaction). As shown in [^, this type of spin- 
spin operator is especially convenient, since it guarantees 
conservation of the quantum number ji in agreement with 
the intrinsic motion of particles in the core possibly being 
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identical to the particles of the three-body system. 

The two-body Hamiltonian is then diagonal in 
the basis {|sfc, jfci)} where ju is the total 

two-body angular momentum after coupling of ji to 
the spin of the core Sk- We denote the correspond¬ 
ing two-body eigenvalues by The natu¬ 

ral basis to describe the three-body system is then 
{\[sk,{siAki)ji]jki,{tki,Sj)jj-,J)}, where J is the total 
three-body angular momentum. The two-body Hamilto¬ 
nian Hki with eigenvalues diagonal in this 

three-body basis. This is in contrast to the remaining 
two-body Hamiltonian Hkj. 

We can calculate the matrix elements of the three-body 
Hamiltonian in this basis. With the expression 

I [Sfc, (si, £ki)ji] jki, i^kj, J) = 

3k j 

? 11 > 4 ) 4 ] ( 5 ., 4 .)jd J) 

Ji ’J Jki J 

I 



we obtain by recoupling that 


A^kididki) 


{qjki\H\q'jki) = 

+ l)(2j;^ + 1) 

X J2(2jkj + 1)4 


( 20 ) 


,Jj ,jkj ) 

"kj 


3kj 


/ jJ 

Sk 

jkj 1 

/ jj 

Sk 

jkj 1 

1 J* 

J 

jki j 

1 

J 

jki J 


where q = {(ki^^kjjjijjj} Therefore the three-body 
Hamiltonian is diagonal in blocks defined by the quantum 
numbers q. 

As an example, we show the 2x2 block correspond¬ 
ing to a core and two light particles all three with spin 
1/2. We further assume ji = jj = 1/2, and J = 1/2 
confining jki and jkj to the values 0 and 1. We note that 
Eq.lHI) depends on the two-body relative orbital angular 
momenta, £ki and £kj, only through the two-body ener¬ 
gies. Then the corresponding block is given by: 


(H) 


2 ^3ki—0) , 1 Tpi^kj T 2 , 3 rpi^kj , 2 dkj—^) 

^ki + i^kj + A^kj 

V 3 ( Tpi^kj , 2 Jfcj — 1 ) Ji^i'^kj ■, 2 dkj — 0 ) \ 

4 k^kj ^kj ) 


, 2 rpi^kj , 2 ’Ifci 

4 Wfcj ^kj ) 

p(^fci,'2 Jfci —1) I 3 Tpi^kj T2 Jkj—^) , 1 ji^i^kj , 2 J kj—^) 
^ki + A^kj + A^kj 


( 21 ) 


This illustrates that even for the schematic case in 
Figd where particles i and j do not interact with each 
other, the three-body Hamiltonian is not separable any¬ 
more. The three-body eigenvalues are then not given by 
the sum of the two-body eigenvalues. When particles i 

and j are identical and £ki=£kj=£, then = 




and the eigenvalues of the matrix ED are 


and + Only the first of these eigenvalues 

corresponds to an antisymmetric eigenfunction under ex¬ 
change of particles i and j. It is precisely twice the av¬ 
erage energy of the two two-body energies. Due to the 
Pauli principle, the two light fermions must occupy both 
the two possible two-body states with angular momen¬ 
tum i. Then the relevant energies are not the individual 
energies of the two two-body states, but their average 
value. 

When none of the two-body interactions contains the 
spin-spin term we have that 


Tpi^ki , 2 dki — 0 ) 
^ki 

j^i^kj , 2 Ofcj —0) 
^kj 


^ki — ^ki 

J^i^kj Jkj—^) _ j^i^kj,^) 
kj kj 


( 22 ) 

(23) 


and the Hamiltonian m is diagonal, with identical di- 

agonal terms given by if/j , which is the 

result obtained in the previous section for spin-zero sys¬ 
tems. When only one of the two-body interactions has 


a non-zero spin-spin term (e.g., Vki, Eq. (1^ is invalid 
and Eq. 123) is valid) the Hamiltonian is diagonal in the 
basis {\ [sk,isi,^ki)ji]jki,{£kj,Sj)jj-,J)} with the ener- 

gies Ef,^ +E^. and +E^j 

There are then two possible three-body states, that again 
are given by the sum of the two-body energies. 


IV. GENERAL PROPERTIES 

The results obtained in the previous section for the 
system in FigUl can be taken as a test for the numerical 
method used to compute three-body states. Taking then 
the schematic system as starting point we analyze the 
main properties of the three-body states, and how they 
evolve in the energy plane when different ingredients are 
added to the calculations. 


A. Bound states and resonances 

We maintain the schematic model where V23=0 and 
mi=oo. We first assume that both V 12 and V 13 only act 
in relative p-waves producing resonances of complex en¬ 
ergies Ei 2 and if 13, respectively. The three-body compu¬ 
tation then leads to the results shown in Fig.|21for a spe¬ 
cific set of parameters. In the small panel, together with 
the discretized complex rotated three-body continuum 
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FIG. 2: Three-body resonance energy and radial wave 
function for 1 ^ 3=0 and two gaussian p-wave interactions 
Vi 2 and Via with ranges 6 i 2 = 1 . 5 fe and 6 i 3 =fe. The 
strengths produce resonances at 2mfe^i?i2/ft^=0.073—i0.033, 
2mb^Ei3/h'^—0.02A—i 0.007. The masses are mi=oo, m 2 =m, 
m 3 = 3 m. Small graph: The solid, long-dashed and short- 
dashed lines connect the computed continuum spectra for 
rotation angles 0=0.23, 0.25, 0.27, respectively. The circle 
above these lines is a three-body resonance independent of 
8. Big graph: The thick lines are the real (solid) and imagi¬ 
nary (dashed) parts of the lowest hyperradial adiabatic wave 
function of the three-body resonance when 0=0.25. The cor¬ 
responding thin (solid, dashed) lines show the asymptotic be¬ 
haviour. 

states (plotted for three different scaling angles), there 
are two additional complex rotated branch cuts starting 
at each of the two-body resonance energies (crosses in 
the figure). These branch cuts correspond to two of the 
particles in a two-body resonant state and the third par¬ 
ticle in the continuum. To keep the figure cleaner these 
cuts have not been plotted. Taking b as an arbitrary 
length unit, the three-body system has a resonance at 
the energy 2mb‘^E/h?=0.097—i0.0A0, clearly distinguish¬ 
able from the continuum background and independent 
of the rotation angle. The computed energy is precisely 
at E=Ei 2 +Ei 3 , as found in section III. This numeri¬ 
cal result is then proving the efficiency of the numerical 
method. To illustrate that this indeed is a three-body 
resonance we also show the lowest radial wave function 
computed for a given rotation angle 9. Both the real and 
imaginary parts vanish asymptotically following the cor¬ 
responding Hankel function 77^+2 required for a pole 
of the S'-matrix as seen from Eq. ©■ 

When one of the two-body systems is bound and the 
other has a resonance the three-body state with energy 
equal to the sum of the two-body energies is simply a 
two-body resonance of one particle relative to the two- 
body bound state. This is illustrated in Fig.|3 where one 
two-body resonance remains unchanged while the other 
two-body state is varied from resonance to bound state. 
When one two-body system is bound, and thus appearing 


2mb^E„/h^ 

0 0.04 R 0.08 0.12 



FIG. 3: Real (Er) and imaginary (Ei) energy {E = Er+iEi) 
of the three-body states (squares) for a system with parame¬ 
ters as in Fig. |21 The resulting energy is the sum of a fixed 
two-body resonance energy (black circle) and the energy of a 
varying two-body state (open circles). When this two-body 
state is bound the three-body state is the two-body resonance 
relative to the bound state. 

on the negative real energy axis (open circles), the three- 
body energies (squares) appear at the same distances to 
the left of the fixed two-body resonance (black circle). 

If the attractive interaction binding the two-body sys¬ 
tem decreases, at some point the two-body state enters 
through the origin into the fourth quadrant of the energy 
plane and becomes a two-body resonance. In parallel, the 
three-body energy (squares) approaches the two-body 
resonance energy (black circle) and continues through the 
two-body resonance in southeastern direction all the time 
following the addition rule. 


B. Virtual states 

The complex scaling method can be viewed as an an¬ 
alytic continuation into complex coordinates, such that 
resonances are “pulled out” of the continuum. Unfortu¬ 
nately the method fails when the pole corresponds to a 
virtual state, due to the necessary large rotation angle. 
Virtual states remain in the unphysical Riemann sheet, 
and numerical investigation of their effects is more diffi¬ 
cult. 

Let now the interactions V 12 and V 13 correspond to a 
virtual s-state and ap-resonance of energies E 12 and E 13 , 
respectively (V23=0, toi=oo). Then, a 1“ three-body S'- 
matrix pole should be present at E=Ei 2 +Ei 3 - To test 
this numerically we include an attractive s-wave interac¬ 
tion V 23 which, combined with the other interactions, is 
sufficiently strong to reveal the existence of a 1 “ three- 
body resonance. We show in Fig. 0] how the complex 
energy of this resonance moves as the strength of V 23 is 
varied. The strongest attraction (the point closest to the 
origin) corresponds to a virtual state in the 2-3 subsys- 


































FIG. 4: Energies of the 1~ three-body states for a system 
with the same masses as in Fig. The p-wave and s-wave 
interactions, Vrs and V12, correspond to a resonance and a 
virtnal state indicated by the crosses. The thick open square 
is the sum of the two two-body energies. Inclusion of a vary¬ 
ing attractive s-wave interaction V23 leads to the three-body 
resonance energies given by the solid circles. For a fixed value 
of V23, when the two-body p-resonance moves as shown by 
the open triangles, then the three-body resonances do it as 
given by the open circles. 


tem at an energy of about —0.004 in FigH A very small 
additional attraction would bind the three-body system, 
which would be Borromean if the three-body system be¬ 
comes bound before the 2-3 subsystem. 

When the strength of V 23 decreases the three-body res¬ 
onance moves (closed circles) towards the energy E 13 of 
the p-resonance. If we could track the three-body reso¬ 
nance for even weaker V 23 we would for some non-zero 
value find that it precisely coincides with E 13 (cross), 
reaching then the discontinuity cut of the Riemann sheet. 
For weaker V 23 the three-body pole could move through 
this cut into any other Riemann sheet, but it is tempting 
to suggest that these states could be interpreted as two- 
body virtual s-states of particle 2 relative to the resonant 
state of particles 1 and 3. Thus such a three-body state 
is not a three-body resonance but a virtual state on the 
unphysical Riemann sheet, although it is not possible to 
see numerically how this energy (square) is reached when 
R23 = 0. 

This can be better understood if the s-state is bound 
as illustrated in Fig.|21 When the s-wave attraction is re¬ 
duced until the bound state energy is zero the three-body 
state moves in parallel towards the resonance energy. A 
continued decrease of the attraction turns the two-body 
bound state into a virtual s-state, and the two-body pole 
moves continuously from the physical to the unphysical 
Riemann sheet. In parallel, the three-body pole moves 
continuously through the p-resonance onto another Rie¬ 
mann sheet. 


The connection between two and three-body reso¬ 
nances are very intimate even when all three masses and 
interactions are non-zero. This is emphasized in FigEI 
where V 13 is varied for fixed ¥ 23 - Then the 1“ three-body 
resonances (open circles) follow precisely the motion of 
the two-body resonances (open triangles). The double 
arrows connect each two-body resonance with its corre¬ 
sponding three-body state. A variation of the two-body 
complex energy produces a similar energy change in the 
three-body state. 

For completeness we also notice that virtual states in 
both subsystems lead to virtual three-body states with an 
energy following the sum rule and sitting in the negative 
energy axis of the unphysical Riemann sheet. 


C. Finite mass and angular momentum coupling 

In a less schematic model also mi has to be finite. Still 
often the three-body structure is mostly influenced by 
two two-body subsystems each dominated by one partial 
wave. Thus we relax the condition mi =00 but maintain 
the reduced masses pi 2 and pi 3 by adjusting m 2 and m 3 . 
The two-body properties then remain unchanged for the 
same interactions V 12 and V 13 . The finite masses destroy 
the separability into two independent subsystems, and 
the eigenvalues of the three-body Hamiltonian are not 
any more given by the sum of the two-body eigenvalues. 
Also, a given combination of partial waves may couple to 
several total angular momenta with different energies. 

When two partial p-waves both contribute, the total 
angular momentum and parity must be 1 “'", 2 ''", 

and combining a p-wave and a d-wave we get L^=l~ , 2 “, 
3“. The corresponding three-fold degeneracy is broken 
for finite mi as illustrated in the upper part of Fig.|3 In 
general, for fixed structure of the system (fixed reduced 
masses, interactions, relative distances between parti¬ 
cles...), the three-body relative kinetic energy is smaller 
for finite mass of the core than for infinite mass, as seen 
by comparing in both cases. Thus 

the three-body resonances often tend to move towards 
the origin. 

However, simultaneous conservation of angular mo¬ 
mentum and center-of-mass might require a change of 
structure resulting in less total interaction energy than 
from the two resonances. Thus, the three-body resonance 
can move in all directions (Fig. [^I. As in the previous 
calculations, only the components corresponding to two- 
body resonances are included. The finite mass effects 
are relatively small in all cases except for two p-waves 
coupled to L'^=l+. This can be understood by express¬ 
ing the Faddeev component related to the 1-3 subsystem 
in the Jacobi coordinates of subsystem 1-2. For L=1 a 
p-wave “rotates” fully into a p-wave. Then the p-wave 
attraction is fully exploited and the resulting three-body 
state has the lowest possible energy. When L and the 
two two-body angular momenta are equal there is a sim¬ 
ilar tendency, decreasing with L, to maximize the total 
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FIG. 5: Upper part: Three-body resonances for a system 
with two-body rednced masses and resonance parameters as in 
Fig-H(p+P case). In the p-\-d case the highest p-resonance has 
been snbstituted by a d-resonance with 2mfe^i5i2/fi^=0.145- 
i 0.009. Lower part: Three-body resonances for a sys¬ 
tem with pi 2 =/ri 3 = 0 . 8 m and p-resonances at Ei 2 =Ei 3 with 
2mfe^i5i2/fi^=0.143-i 0.059. The fixed two-body resonances 
are indicated by the crosses and the three-body resonances by 
a big open circle (mi=cx 3 ), and solid circles (L= 0 ), squares 
(L=l), triangles (L=2), and stars {L—3) for finite mi values. 


attraction. Other components obtained by rotation are 
less important and uneven two-body angular momenta 
cannot exploit the attraction as efficiently, producing a 
much smaller change with mi. In practice detailed quan¬ 
tum mechanical computations are needed for each case 
to determine the size of the effects. 

It is remarkable that the effects of the coupling and the 
center-of-mass motion can lower the energy sufficiently 
to produce a 1+ three-body bound state (lower part of 
Fig.EJ). This is a Borromean system obtained by only two 
two-body interactions. When one resonance is replaced 
by a virtual state we have not numerically been able to 
find any three-body resonance produced entirely by finite 
mass effects and without change of interactions. 



FIG. 6 : Evolution of the three-body resonances for ®He (cir¬ 
cles), ®Li (squares), and ®Be (triangles), when the nucleon- 
nucleon interaction is progressively introduced. The big open 
symbols correspond to zero nucleon-nucleon interaction. The 
0 ^ and 2 "*" resonances are shown by the closed and open sym¬ 
bols respectively. The stars in the ®Li case are the I’*' res¬ 
onance corresponding to the T = 0 channel in the neutron- 
proton interaction. The crosses show the alpha-nucleon reso¬ 
nance energies. 


V. NUMERICAL ILLUSTRATIONS 

The origin of the resonances and bound states can 
be illustrated by different examples connecting the bare 
schematic model with well known realistic nuclear struc¬ 
tures. We first study two nucleons outside the spin-zero 
a-particle core, then we present detailed realistic compu¬ 
tations for the same A = 6 systems, and finally we extend 
to the non-zero core-spin of the Borromean nucleus ^^Ne. 


A. Core with zero spin 

Let us consider a system made by an infinitely heavy 
core and two non-interacting spin 1/2 particles. The core 
is assumed to be charged (with twice the proton charge) 
and the particle-core reduced mass is 0.8m, where m 
is the nucleon mass. This system is very similar to an 
Qf-particle and two nucleons. The particle-core reduced 
mass is the same, and the center-of-mass effects are very 
small and not relevant for our present purpose |3l|. 
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Let us also consider a simple particle-core interaction 
given by a central short-range p-wave interaction pro¬ 
ducing a particle-core (uncharged) resonance at an en¬ 
ergy of 0.74 MeV with a width of 0.60 MeV. When the 
light particle is assumed to have the proton charge, the 
same short-range interaction produces a core-particle res¬ 
onance with energy and width of 1.60 MeV and 1.26 
MeV, respectively. These values are consistent with the 
lowest p-resonances in ®He and ®Li, respectively 

In Fig El we show the results for one heavy and two 
light particles resembling the ®He, ®Li, and ®Be systems. 
The upper, middle and lower parts correspond to both 
uncharged, one neutral and one with the proton charge, 
and both with the proton charge, respectively. When the 
nucleon-nucleon interaction is equal to zero, the three- 
body system must have a resonance at an energy equal 
to the sum of the two two-body resonance energies. This 
is again confirmed by calculations using the complex ro¬ 
tated adiabatic expansion method. We obtain the three- 
body resonances indicated by the big open circle, big 
open square, and big open triangle in the upper, central, 
and lower parts of Fig El They match precisely with the 
sum of the two-body energies each indicated by a cross. 

The known spectrum for the schematic three-body sys¬ 
tem, determined completely by the internal two-body 
states, can be used as starting point to trace the connec¬ 
tion with the “realistic” three-body system. This must be 
done by including the effects produced by i) an additional 
interaction between the two light particles, ii) quantum 
numbers conservation (like the angular momentum), and 
iii) center-of-mass effects originating from the finite mass 
of the core. 

Effects i) and ii) are closely connected. Once the in¬ 
teraction between the two light particles is included the 
three-body wave function, initially independent of the 
total angular momentum (Eg. 11211 1. is now depending on 
L. This can be seen in the upper part of EigEl where 
we include the neutron-neutron interaction multiplied by 
a global factor varying from 0 to 1. Then the 0+ and 
2 + states evolve as shown by the solid and open circles, 
respectively. When the full neutron-neutron interaction 
is included (last point on each curve), the system similar 
to ®He has a bound Borromean 0+ state (with a bind¬ 
ing energy close to —1.9 MeV) and a very narrow 2+ 
three-body resonance with energy 0.34 MeV. 

For the system similar to ®Li shown in the middle part 
of the figure one of the light particles has the charge of 
the proton. The neutron-proton interaction is again con¬ 
tinuously switched on from zero to full strength. The 
0+ and 2+ states (T=l), analogous to those in the up¬ 
per part, move as given by the closed and open squares, 
respectively. The final 0+ state is still below the three- 
body threshold, with a binding energy of —1.0 MeV. The 
final 2+ resonance has an energy of 1.1 MeV and a width 
of 0.1 MeV. The shift in energy compared to the states in 
the upper part is produced by the Coulomb repulsion be¬ 
tween the core and the charged particle. Now the allowed 
T=0 coupling gives rise to additional three-body states. 



Scheme 

Real, (no 3-b) 

Realistic 

Exper. 

®He 

0 + 

-1.87 

-0.04 

-0.96 

-0.97 ±0.04 


2 + 

(0.34,0.01) 

(1.36,1.12) 

(0.87,0.11) 

(0.83,0.11) 


0 + 

-0.93 

(0.75,0.15) 

-0.14 

-0.14 

®Li 

1 + 

-5.60 

-3.05 

-3.73 

-3.70 


2 + 

(1.07,0.09) 

(1.92,0.87) 

(1.66,0.50) 

(1.67,0.54) 

®Be 

0 + 

(0.55,0.16) 

(2.05,0.58) 

(1.37,0.11) 

(1.37,0.09) 


2 + 

(2.32,0.46) 

(3.10,1.89) 

(3.02,1.65) 

(3.04,1.16) 


TABLE I: Computed and experimental energies (in MeV) for 
®He, ®Li, and ®Be. The energies within brackets are {Er, Bij), 
that give the resonance energy and its width. The computa¬ 
tions correspond to the schematic system described in the 
text (second column), realistic calculations without inclusion 
of an effective three-body force (third column), and the same 
calculation when a three-body force is used (fourth column). 
The experimental values (last column) are taken from 1321 . 
When not specified, the error bars are smaller than the last 
digit. 


e.g. the I"*" state shown by the stars. When the three- 
body 1+ state becomes bound the system is still Bor¬ 
romean. However, for some threshold strength the two 
light particles become bound, and the system is not Bor¬ 
romean anymore. When the full interaction is included 
the two light particles then form a deuteron nucleus, and 
the three-body I'*' state becomes the ground state of the 
system, very similar to ®Li, with a binding of about —5.5 
MeV (outside the scale of the figure). 

In the lower part of the figure we show results for two 
light particles each with the proton charge. As in the up¬ 
per part only the T = 1 states are allowed, and since all 
the three two-body subsystems feel the Coulomb repul¬ 
sion, the energies of the 0+ states (closed triangles) and 
the 2+ states (open triangles) are clearly larger than in 
the previous cases. In fact, when the full proton-proton 
interaction is included the O'*" ground state is unbound, 
with resonance energy and width of 0.5 MeV and 0.2 
MeV, respectively. For the 2+ resonance the energy and 
width are 2.3 MeV and 0.5 MeV. 

The effects produced by the finite mass of the core 
(center-of-mass effect), do not change the previous ener¬ 
gies significantly (typically no more than 100 keV [^1. 
The spectra obtained for the “simplified” ®He, ®Li, and 
®Be nuclei (only a simple p-wave particle-core interaction 
has been used), indicate where to find the true states. 
This is seen in tabled where the second column gives the 
computed energies. A comparison with the experimental 
results (last column) reveals that in all the three nuclei 
the main features of the spectrum are well reproduced. 
Therefore, starting from the schematic case where the 
three-body energy is given by the sum of the two-body 
energies and taking into account the main characteristics 
of the remaining two-body interaction, it is possible to 
estimate where the three-body states must be placed. 

As seen in the table, the computed energies are always 
lower than the experimental ones. This is because the 
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nucleon-core interaction is without any spin-orbit term, 
and consequently the two possible p-resonances (pi /2 and 
P 3 / 2 ) appear at the same energy. This is overbinding the 
system, since the pi/ 2 -i'esonance is known to be a few 
MeV higher , and a realistic calculation must include 
a repulsive term pushing up this resonance to the cor¬ 
rect value. On top of this, other than p-waves should 
be included in the calculation. For s-waves the effect 
of the Pauli principle must be accounted for, since the 
Si/ 2 -shell is fully occupied in the a-core. This effectively 
amounts to a highly repulsive s-wave potential at short 
distances, also pushing the computed energies towards 
higher values. Furthermore, it is well known that three- 
body calculations using only realistic two-body interac¬ 
tions typically underbind the system. This can be cured 
by inclusion of an additional potential taking into ac¬ 
count possible three-body effects to reproduce the exper¬ 
imental values. 

Realistic detailed calculations using the (complex 
scaled) hyperspheric adiabatic expansion method con¬ 
cerning ®He can be found in plilL l34l|. For completeness 
we show in the next section the results obtained for ®Li 
and ®Be when similar calculations are performed. 


B. Realistic calculations for ®Li and ®Be 

In the use of phase equivalent potentials is sug¬ 
gested as an efficient method to take into account 
the Pauli principle in three-body calculations. This 
method, together with the hyperspheric adiabatic expan¬ 
sion method, is used to compute the ground state of ^^Li 
and ®He. Thus, for ®He the alpha-neutron s-wave inter¬ 
action is able to bind the neutron into a Pauli forbidden 
state. To exclude this state we use a phase equivalent 
potential with exactly the same phase shifts for all ener¬ 
gies, but with one less bound state. The p-wave interac¬ 
tion contains a central and a spin-orbit term, such that 
the P 3/2 and pi /2 resonances in ®He are placed at the 
experimental energy values |32|. The neutron-neutron 
interaction can also be found in |^. 

The calculation using only the two-body interactions 
underbinds the three-body nucleus (third column in ta¬ 
ble ni, and an additional effective three-body force is 
needed to recover the experimental value (fourth col¬ 
umn). In the same method supplemented by com¬ 
plex scaling is used to compute the 2+ resonance in ®He. 
In the present work the calculation is slightly simplified 
by using a repulsive s-wave alpha-neutron interaction to 
take the Pauli principle into account. When the same 
phase equivalent potential as for the ground state is used, 
the computed 2+ state in ®He is also given in the third 
and fourth columns of tabled As before, the calculation 
without a three-body force is underbinding the system 
by about 0.5 MeV, and an effective three-body force is 
again needed. In all the calculations the three-body force 
is assumed to be gaussian with a range of 3 fm and the 
strength is adjusted to fit the experimental value. For 


bound states it is obviously possible to find a strength 
fitting precisely the experimental value. For resonances, 
this single parameter is fitting simultaneously the energy 
and the width of the resonance. 

For ®Li and ®Be we perform exactly the same calcu¬ 
lations but adding the corresponding Coulomb interac¬ 
tion for the two cases. The additional proton-core repul¬ 
sion in ®Li is pushing up the energies of the 0+ and 2+ 
states compared to the energies obtained for ®He, and 
even more for ®Be, where all the three two-body subsys¬ 
tems feel the Coulomb repulsion. The computed energies 
are again given in the third column of tabled when only 
the two-body interactions are used. As before, an effec¬ 
tive three-body interaction is needed to fit the experi¬ 
mental energies (fourth column). The O'*' state in ®Li is 
still below the two-nucleon separation energy threshold, 
although when the three-body interaction is suppressed 
the computed state appears to be unbound. Again, when 
the strength of the three-body force is used to match the 
experimental resonance energies, the widths are also well 
reproduced, except for the 2+ state in ®Be, where the 
width is significantly larger than the experimental value. 

In addition to the O'*' and 2+ T = 1 states, in ®Li it 
is also possible to have states with T = 0. These three- 
body states correspond to structures where the neutron 
and the proton form a deuteron (d) nucleus, suggest¬ 
ing a description of the nucleus as a d+a two-body sys¬ 
tem. Therefore, we have performed a three-body calcu¬ 
lation, identical to the ones described above, but where 
the nucleon-nucleon interaction has been substituted by 
the T = 0 interaction used in |^. This potential re¬ 
produces the experimental deuteron binding energy, root 
mean square radius, electric quadrupole moment, and 
provides a d-wave content of 5.6%. The correspondig 
computed energies for the 1 + (ground) state are also 
given in tabled 

The three-body Thomas-Ehrman shifts of these 
isobaric analog O'*" and 2 + states are then obtained from 
tableE i.e. in MeV (1.40,1.17,1.05) and (1.83,1.80,1.65) 
for the realistic cases without and with three-body poten¬ 
tial, respectively. The latter results coincide by definition 
with the experimental values. The decreasing tendency 
in the computed differences are due to the smaller en¬ 
ergy and the larger radii which in turn leads to smaller 
effects of the Coulomb interaction. The underbinding is 
responsible for this discrepancy with measurements. 


C. Core with non-zero spin 

When the three particles have non-zero spins and spin- 
dependent two-body interactions, as in any realistic nu¬ 
clear potential, then the three-body Hamiltonian is not 
even separable in the schematic model with an infinitely 
heavy core and two mutually non-interacting particles, 
see section null The Hamiltonian matrix is then or¬ 
ganized in blocks with the matrix elements given by 
Eg.CT. 
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FIG. 7: energies for the 1/2”, ^ and 5/2” states as 

a function of the percentage of the proton-proton interaction 
included in the calculation. 

An example is for ^^Ne, that is well described as a 
three-body system made by an core and two pro¬ 
tons 1^. The core of the system has negative parity 
and spin 1/2. The two-body subsystem, (^®0-|-p), 
has four low-lying resonances, two of them ( 0 ” and 1 ”) 
arise from the coupling of a relative Si/ 2 -wave and the 
spin of the core, and other two (2” and 3”) result from 
the coupling between a (i 5 / 2 -wave and the spin of the 
core. The experimental energies and widths of 

these two-body resonances are (0.535,0.040±0.020) MeV, 
(0.728±0.006,(0.040) MeV, (0.959±0.005, 0.040±0.030) 
MeV, and (1.256±0.004, (0.004) MeV, respectively 

In several proton-core interactions, all of them 
reproducing the experimental energies and widths, are 
given. These interactions contain spin-dependent opera¬ 
tors, in particular the spin-splitting operator s^-jp, where 
Sc is the spin of the core and jp results from the coupling 
between the relative orbital angular momentum and the 
spin 1/2 of the proton. In the following calculations we 
shall use the s and d-wave interactions corresponding to 
a gaussian two-body potential with the parameters given 
in Table 1 of [s^l, but where the spin-orbit strength has 
been modihed to push away the d 3 / 2 -resonances in 
but keeping the two (i 5 / 2 -i'esonances at the right energy. 
In this way our slightly simplihed ^^Ne is characterized 
exclusively by the 0”, 1”, 2”, and 3” resonances in ^®F, 
and we can perform for ^^Ne the same kind of analysis 
as in the previous cases. 

Let us begin by performing calculations for the 
schematic model without any proton-proton interaction 
and for an infinitely heavy core. The effect of the finite 
mass of the core is clearly smaller than for ®He, ®Li, and 
®Be, and the core can safely be assumed to be infinitely 
heavy. 

When both protons are in an s-wave relative to the 
core, only the 0 ” and 1 ” states in ^®F are obtained. 


and only the components with relative core-proton or¬ 
bital angular momentum equal to zero are included in 
the calculation. In this case the total three- body angular 
momentum can only take the value J = 1/2 |37| . and the 
block to be diagonalized is given in Eq. mi- From the two 
possible eigenfunctions only one of them is antisymmet¬ 
ric under exchange of the two protons, and corresponds 
to the eigenvalue + 3E^^ ^)/ 2 , that is twice the 

average energy of the two Si/ 2 -resonances. This gives rise 
to a 1/2” energy of 1.32 MeV, that is recovered when the 
calculation using the complex scaled hyperspheric adia¬ 
batic expansion method is performed. 

When the proton-proton interaction Vpp is progres¬ 
sively introduced, the contribution from the d-waves 
starts to be relevant, and the components with rela¬ 
tive core-proton orbital angular momentum equal to 2 
are also included. The energy of the 1/2” state then 
changes as shown by the solid circles in FigCI Around 
90% of the proton-proton interaction is already binding 
the three-body system, becoming then Borromean. The 
full proton-proton interaction gives a binding of about 
0.40 MeV and a d-wave content of 37%. This binding is 
0.3 MeV smaller than the one given in table 8 of |3^ . 
This is due to the fact that in this calculation, in order 
to clean the structure of ^^Ne, we have artificially pushed 
up the d 3 / 2 -resonances by using a large spin-orbit split¬ 
ting. This is also reducing the d-wave content compared 
to the result in |^. 

In ref.|^ we show how the 3/2” state in ^^Ne is dom¬ 
inated by the sd-interferences between the components. 
To investigate this resonances both (0 and 2) core-proton 
relative orbital angular momenta are then needed and 
therefore included in the calculation. For the schematic 
^’^Ne model (Vpp = 0), one of the protons can be in one 
of the s-resonances in ^®F (0” or 1”), and the other can 
be in one of the d 5 / 2 -i'esonances (2” or 3”). Assuming 
J = 3/2, the block in Ea. lEHIl is already diagonal with 
only one non-zero energy equal to E^^ '> -|- E^"^ \ This 
energy of 1.65 MeV is recovered in the numerical cal¬ 
culation, which also reveals that the only non-vanishing 
components correspond to one proton in the s-wave I” 
resonance and the other in the d-wave 2” resonance. The 
evolution of the energy when the proton-proton interac¬ 
tion is introduced is shown by the white circles in FigEI 
reaching a hnal value of 0.94 MeV. This value is also 
around 0.3 MeV less bound than the one given in table 
8 of [ 33 . 

Finally, in the figure we also show the evolution of 
the lowest 5/2” resonance in ^^Ne (solid squares) when 
Vpp is introduced. This state is also dominated by sd- 
interference components. 

VI. QUALITATIVE CONSIDERATIONS 

Using the schematic case in which V 23 = 0 as starting 
point, it is possible to trace the behaviour of the three- 
body resonances (Figs. [3 El andEJ when the center-of- 
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mass effects and the V 23 interaction are included. Taking 
into account the main features of the 1^23 interaction, 
it is then possible to make crude estimates concerning 
different realistic systems. 

For nuclear systems we take the mass unit m equal 
to the nucleon mass and a length unit equal to a typi¬ 
cal range for the nuclear interaction (6 = 2 fm). With 
this choice, the lowest p-resonance in the neutron-alpha 
interaction has dimensionless energy 2mb‘^\E[i\/f{} « 
0.15. Therefore, in ®He (a-|-n-|-n), before including the 
neutron-neutron interaction, the three-body system must 
have a resonance with dimensionless energy of about 0.30. 
The neutron-neutron interaction has a low lying virtual 
state at 2mh'^\ER\/ti^ « 0.02, which is enough to bind 
the three-body state as evidenced by the bound state in 
®He at 2mb^Eii/h^ « —0.19. 

For ^^Li (®Li-|-n-|-u) the s-waves are pronounced and 
the lowest neutron-®Li virtual s-state would after spin 
averaging appear at about 0.07 which, combined with 
the neutron-neutron interaction, is sufficient to bind the 
three-body system. Substituting ®Li by a A-particle we 
arrive at an unstable system. With the A-nucleon s- 
wave scatterin g le ngth of 2 fm obtained from hypertriton 
computations the virtual s-state would be at around 
1.7, too large to bind the A-neutron-neutron system. 

We can continue to a more speculative system obtained 
by combining the recently discovered pentaquark reso¬ 
nance with a third particle like a meson. The pentaquark, 
9~^ ^ has a mass of 1540 MeV and decays mainly into a 
kaon and a neutron . The energy above threshold 
is then around 107 MeV, that corresponds to a dimen¬ 
sionless energy of about 20. Thus, the neutron-neutron 
interaction, that was able to bind a-l-n-l-n and ^hi+n+n 
respectively with 0.15 and 0.07, but not A-|-n-|-n with 
1.7, cannot bind the K'^ + n + n system. Thus a nucleon- 
kaon-meson cluster system would also exploit the in¬ 
trinsic quark degrees of freedom and not resemble any 
three-body structure although additional binding would 
be picked up. In any case the narrow width of 9'^ strongly 
suggests that the resonance cannot have nucleon-kaon 
(two-particle) structure, but rather is a genuine five- 
quark system or perhaps a more exotic structure. 

The semi-quantitative knowledge obtained with the 
schematic model can also be used for molecules due to 
the use of dimensionless parameters. For molecules we 
choose as typical scale units the mass of the ^He atom 
and the range 6 = 10 A. The ^He-l-'^He dimer has then 
2mb'^\Eji\/h'^ « 0.75 We can then expect that, ex¬ 
cept for center-of-mass effects, the trimers ^He-|-'*He-|-"'^He 
and ^He-|-^He-|-^He should both have an unbound state 
with dimensionless energy of about 1.50 when the inter¬ 
action between the two identical particles is neglected. 
However, including these interactions of dimensionless 
energies 0.02 for ^He-|-"'^He and 5.2 for ^He-|-^He, we ar¬ 
rive at trimers ^He-|-'*He-|-"’^He and ^He-|-^He-|-^He respec¬ 
tively bound and unbound as observed. 


VII. SUMMARY AND CONCLUSIONS 

The two-body interactions are assumed to determine 
completely the three-body structure including the con¬ 
tinuum properties. Close to the threshold of binding, 
slightly above or below, the structure of the three-body 
system mostly depends on the two-body low-energy scat¬ 
tering properties. These properties are periodically re¬ 
peated with increasing strengths of these attractive po¬ 
tentials corresponding to one or many bound states in 
each of the investigated channels. Therefore it is suffi¬ 
cient to study unbound or weakly bound two- and three- 
body systems. 

The low-energy scattering properties are determined 
by the phase shifts and reflected in the poles of the 
5'-matrix, i.e. by the resonances and virtual states. 
Thus, changing the two-body interactions simultaneously 
change the energies of both two- and three-body reso¬ 
nances. The relative changes of these energy observables 
are expected to be intuitively easier to understand than 
by using the connection via the two-body interaction 
strengths. The strategy is then to vary the two-body 
resonances and virtual states and study the changes of 
the corresponding quantities in the three-body system. 

We employ the efficient and well tested method of 
complex scaled hyperspherical adiabatic expansion. This 
method is first briefly sketched to define the notation. 
Then we define a schematic model with an infinitely 
heavy core and two mutually non-interacting particles. 
We prove mathematically for spin-zero particles that our 
formulation provides a pole in the three-body S'-matrix 
if and only if the complex energy is equal to the sum 
of two complex energies each corresponding to poles of 
different two-body S'-matrices for the two particle-core 
subsystems. The generalization to non-zero spins are for¬ 
mulated and shown to involve diagonalization of simple 
block-diagonal Hamiltonians. 

The general properties are demonstrated numerically 
in schematic examples involving both resonances and vir¬ 
tual states. We then investigate the sizes and trends 
resulting from a finite core mass, non-zero interaction 
between the light particles and the coupling of different 
orbital angular momenta. We use dimensionless units 
to allow easy application on physics systems of differ¬ 
ent scales. Each of these effects can be substantial but 
the structure of the states can be uniquely traced back 
to the origin in the structure of the schematic model. 
Borromean systems can arise with only two interactions 
from center-of-mass effects and a favorable coupling of 
two angular momenta. 

For further numerical illustration we use essentially re¬ 
alistic nuclear three-body systems, ®He, ®Li, and ®Be, 
consisting of two nucleons and an a-particle to trace back 
their measured spectra to our bare schematic model. The 
effects of isospin symmetry and the mixing short-range 
and Coulomb interactions are then seen. For complete¬ 
ness we also present unpublished and fully realistic calcu¬ 
lations of ®Li and ®Be. The origin of the structure is still 
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apparent, but to get accurate energies and wave func¬ 
tions we must include effects of spin-orbit couplings and 
the Pauli principle. Three-body Thomas-Ehrman shifts 
can then be studied for these isobaric analog states. Fi¬ 
nally the effects of non-zero core spin is investigated for 
the Borromean nucleus ^^Ne consisting of and two 
protons. 

In the last section we explain, and illustrate by exam¬ 
ples, how to make qualitative estimates of the three-body 
energies and their structure from the two-body proper¬ 
ties of the subsystems. First we test by known examples 
of Borromean halo nuclei, hypertriton, molecular helium 
clusters and the recently highlighted more speculative 
pentaquark. 

In conclusion, we have demonstrated the strong cor¬ 
relation between two- and three-body resonances. The 
three-body energy and structure can be traced back to 
the properties of an infinitely heavy core and two non¬ 
interacting light particles. Substantial changes are often 
needed to arrive at accurate and realistic properties but 
the generic origin is apparent and revealing both struc¬ 
tures and energies. 


APPENDIX A: USEFUL INTEGRALS 

The three-body ^'-matrix for the schematic model in 
Ea. (ll7ll is obtained through the definition in Ea. IlIfil) by 
inserting Eq. m into Eq. Then the following inte¬ 


grals are needed: 

^7r/2 ^ ^ 

/ da(sin (cos ^ ^ (cos 2a) 

Jo 

xh''l\ka:x)h‘^ljikyy) = (Al) 

(-I)"'7r(sina„)^*(cos 

[kpY 


f dQ;(sin (cos 2 ) 2a) 

Jo 

xhYJ (,kxx)hY\kyy) = 0 (A2) 


f dQ;(sin (cos ^ ^ (cos 2a) 

Jo 

xhYJ(k^x)hYJ(kyy) = 
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da(sm a)^*+^(cos ^ ^ (cos 2a) 
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